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We present a new approach to coarse-graining of variables describing dust flow in GR. It is 
based on assigning quasi-local shear, twist and expansion to 2-dimensional surfaces with the help 
of isometric embeddings into the 3-dimensional Euclidean space and deriving the time evolution 
equations for them. In contrast to the popular Buchert's scheme it allows to coarse-grain tensorial 
quantities in a coordinate-independent way. The framework can be used to estimate backreaction 
in inhomogeneous cosmological models. 
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I. INTRODUCTION 



o 

In this paper we discuss an alternative to the well-known Buchert's averaging scheme [l| for inhomogeneous cosmo- 

■ logical models with dust. The new procedure allows for coarse-graining of both scalar (expansion) and tensorial (shear 
and vorticity) part of the fluid velocity gradient in a coordinate-independent way. Given a 3-f 1 decomposition of the 
spacetime we assign the coarse-grained expansion, shear and vorticity to three-dimensional comoving domains of the 

O fluid. The resulting non-local quantities depend on the 3-1-1 splitting, but not on the coordinate system introduced 
on the spatial slices. If the fluid is irrotational, there exists a preferred, orthogonal splitting and the formalism is 
^ [ entirely covariant; otherwise one needs to fix the constant time slices by an additional condition. 

The procedure is based on a reformulation of the definition of the volume average of the velocity gradient over a 
comoving domain in Newtonian theory. It turns out that it is possible to reexpress that volume average as a functional 
^-H of the time derivative of the metric tensor induced on the domain's two-dimensional boundary. It has been noted 
J> ' before 0, Q that such volume-averaged velocity gradients satisfy an evolution equation which is very similar to the 
I evolution equation of the local velocity gradient, the only difference being a surface integral of the inhomogeneities of 
: the flow over the domain's boundary. This volume average can be rewritten as a surface integral of the fluid velocity 
due to the divergence theorem. The key observation is that the surface values of the velocity field can be almost 
uniquely reconstructed from the comoving time derivative of the metric induced on the domain's boundary as long 
as the that boundary has the topology of S*^ and its metric has a positive curvature. The only part of the velocity 
field which cannot be reconstructed is the total linear and angular velocity, responsible for translations and rotations 
of the surface as a whole. However, since they do not contribute to the average shear and expansion, these quantities 
can be evaluated just from the reconstructed velocity field. Thus it is possible to reexpress in Newtonian theory the 
average shear and expansion as functionals of the derivative of the induced metric on the boundary. 

The reformulated definition carries over easily to the relativistic case. If we fix a comoving, three-dimensional 

■ domain in a 3-f 1 decomposed spacetime, we can find the isometric embeddings of the domain's two-dimensional 
' boundaries into the three-dimensional Euclidean space and evaluate the coarse-grained expansion and shear. If the 

background metric is not fiat, the resulting quantities are not equivalent to any volume integrals any more, but 
nevertheless they can serve as a measure the deformation of a finite fluid element. As we show in the paper, the new 
definition satisfies simple reasonability conditions: it yields the right answer in the limit of the domain shrinking to 
a point and in FLRW metrics. 

We also present a way to define the coarse-grained vorticity from the projection of particle four-velocities to the 
spatial constant-time slice if the dust is rotating. The definition is independent of the coordinates introduced on the 
spatial slice, but in contrast to the irrotational case we do not have a simple, unique definition of a preferred 3-1-1 
splitting and the coarse-grained quantities do depend on it. 

Finally we present the evolution equation for the coarse-grained velocity gradient and show that its structure is 
very similar to the evolution equation for its local counterpart. Two additional terms in the form of surface integrals 
may serve as a measure of backreaction. i.e. the impact of the inhomogeneities of the flow on the large scale motion 
of the dust within the domain. 

The proposed procedure is more mathematically involved then Buchert's, as it requires finding the isometric em- 
beddings of S'^ surfaces, which is equivalent to solving a system of coupled, non-linear PDE's. In most cases this 
is impossible to do analytically, but numerical methods have already been developed for that problem 0, Q • Let us 
also note that in contrast to Buchert's approach we do not have evolution equations for the coarse-grained curvature 
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and matter density. The dust was chosen as the matter model in this paper just for its simphcity; the formahsm 
introduced here can be apphed to any other fluid in General Relativity with only minor modifications. For a more 
exhaustive treatment of the topic, including rigorous derivations of the results, see 



II. COARSE-GRAINING VIA ISOMETRIC EMBEDDINGS 



Consider gravitating, Newtonian, presureless dust in a flat space. It is described by the velocity field u*, density 
p and the Newtonian potential (p. The gradient of the velocity field Q'^ = ^ describes how an infinitesimal fluid 
or dust element is deformed and how it rotates during its motion. It can be decomposed into the scalar expansion, 
symmetric traceless shear and antisymmetric vorticity 

Qij = 2 ^ % + '^{ij) + ^[ij] ■ 
Qij satisfies the evolution equation along the fiow 

^Q^, = -Q\Q',~r,, (1) 

where ^ is the convective derivative (see @, [^@). This equation has a coarse-grained counterpart which one obtains 
by taking the volume average of Q^j over a comoving domain Gt 

Note that (Q* ) is effectively a surface integral over the boundary of dGt 



(3) 



Its evolution equation has a similar structure to ([T]), 

|(Q^) = -(0\)(Q^)-(<^'^> + i?^ (4) 

with (0'* j) being the volume average of the Hessian of the potential over Gt- The additional term 



= 1 / {Sv''6v%nk-Sv''Sv\n,)da, (5) 



measures how the inhomogeneities of the fluid velocity field dv"^ = — S'ff^ct its large scale motion in the 

domain Gt 

In the relativistic case, on a non-flat background, we do not have a well-defined, coordinate-independent notion 
of volume averages of tensor objects, so we cannot use ([2]) to obtain a relativistic counterpart of Q. It turns out, 
however, that definition ([2|) can be reformulated in a geometric, coordinate independent manner which has a natural 
extension to non-fiat backgrounds. We will first deal with the symmetric part of (Qij), i.e. shear and expansion. We 
will show that under certain technical assumptions they can be reexpressed as functionals of the metric induced on 
the boundary dGt and its convective (comoving with the dust) time derivative. 

The main tool we are going to employ is the isometric embedding theorem for surfaces of topology, conjectured 
by Weyl and proved by Lewy, Alexandrov, Pogorelov, Nirenberg and Cohn-Vossen (see 0, Hi): 

Theorem II. 1 (Isometric embedding theorem for S^) Given a compact, orientable surface S homeomorphic to 
S'^ , with positive metric q whose scalar curvature R > 0. Then 

• there exists an isometric embedding 



into the 3-dimensional Euclidean space; 



3 



• the embedding is unique up to rigid rotations, translations and reflexions. 

If we fix the orientation, the theorem states that for a given S*^ surface with a positive metric of positive curvature 
can be recognized as a submanifold in E'^ and this can be done only in one way, up to moving the surface around and 
rotating it as a whole. Surfaces satisfying the hypothesis of theorem lII.il will be called admissible. 

Consider now the boundary dGt of the coarse-graining domain, parametrized by two coordinates 6"^ labeling 
individual particles. It is described, together with its time evolution, by equations 



The induced metric has the form of 



qMt,0^)^C,AC\BS.3 (6) 



and its time derivative in coordinates 6^ 



qAB = 2w\ 4 B^ 



Assume now that dGt is admissible at least for some time and that we know its position in space at one instance 
of one time t = tg. We will show that in that case {Q(ij)) can be read out solely from the induced metric even if 
we do not have each individual particle's velocity and position. 

Given qAsit) we can find the isometric embeddings of dGt, consistent with its orientation, into E'^at each time t. 
We assume that at to the embedding is known, but at other times we need to solve the following system of PDE's 

qAB ^x\ax\b^v C^) 
for three functions x'(i, 9^). We construct this way a one-parameter family of embeddings 

x^^x'(t,0^) 

for which we assume the additional condition 

The first statement of theorem III. II assures that functions x'(t, 0^) can always be found, although they are non-unique. 
The true particle positions of dGt, given by C,'^{t,9^), constitute trivially a family of isometric embeddings, so from 
the the uniqueness part of III. II we know that functions and C a-re related via 

xHt,o^) = R',it)c\t,0^) + w\t), 

where i?*^ is a rotation. If we take a time derivative of this equation at < = <0i we obtain 



dx' 
dt 



with constants $7^ = —^ji and . Thus the velocity field u\ reconstructed from an arbitrary family isometric 
embeddings, differs from the true one only by the restrictions to dGt of a rotational (51*^- x^) and s translational (i?*) 
vector field. This difference does not infiuence the value of integral ^ symmetrized in the two indices 

U(rnj)Aa^l do- = (8) 

dGt JdGt 

so we can use the reconstructed u instead of v to evaluate the shear and expansion. 

In fact we do not have to find the embeddings explicitly to obtain the reconstructed velocity field. By differentiating 
([7]) with respect to t we obtain an equation relating directly to the convective time derivative of qAB at i = to 

'^^\,A C\b)^^3 =qAB- 

The linear operator V(^[Y]ab ~ 2y| ^ 6ij appearing on the left hand side yields the variation of the surfaces' 
induced metric if the embedding itself is being dragged by a vector field . By the virtue of theorem III. II for an 
admissible surface it has a six-dimensional kernel consisting of rigid motions and a "non-unique inverse" V^^. We 
can now rewrite the definition of the symmetric part coarse-grained velocity gradient explicitly as a functional of qAB 

(Q{v))^T7 I ^rM9As](, »^j)dcr. (9) 

^ JdGt 



FIG. 1: Coarse-graining of Z^^ over the slices of a tube of particles 



III. VELOCITY GRADIENT IN GENERAL RELATIVITY 

Formula ^ is clearly more complicated that ([2]), but its main advantage lies in the fact that it can be generalized 
to fluids on a non-flat background in a simple and natural way. Consider a (possibly curved) spacetime, filled with 
dust described by density p and four-velocity field u. The integral curves of u are geodesies, i.e. V^u'^ = 0. The 
four-velocity gradient Z^^^i, = V^Uf^ describes the deformations of an infinitesimal fluid element in its motion. It 
is orthogonal to u in both indices and therefore effectively a three-dimensional object satisfies the relativistic 
counterpart of H]) 

V„Z^, = -Z^^ Z"^ R\^p u^. (10) 

Unless stated otherwise, in this paper we will assume the flow to be irrotational, i.e. Z^i^ = Z{fii^)- We will now show 
how to derive the coarse-grained counterpart of (|10p in which (Zij) corresponds to a finite fluid volume rather than 
an infinitesimal element. 

In comoving and orthogonal coordinates {t,y^) one can perform the standard ADM decomposition of the metric 
into the spatial part hij, with vanishing iV* and = 1 due to the coordinate system choice. Consider a three- 
dimensional tube of particle worldlines dC enclosing a cylinder C (see fig. [T|). The tube is described by equations 

= C being three time-independent functions and 6^ comoving coordinates on the tube labeling individual 

particles. The intersections of the cylinder with constant time slices St, denoted in this paper by dCt, constitute two- 
dimensional surfaces of spherical topology. Assuming that the surfaces are admissible we assign the coarse-grained 
velocity gradient to their interior in the following way: we take the induced metric 

qAB{t,9^)^C.Ae,Bh^J{t) (11) 

on dCt and find an arbitrary family of isometric embeddings ft : dCt dDt C E'^. They are described by equations 

a;" = x"(i,^^) 

with three functions satisfying a system of non-linear PDE's 

gAB(t,e-^)-X%X%'5ab. 

These are subject to gauge transformations of the form of rigid motions 

x'^it, e^) ^ R\{t) x\t, e^) + w-{t). (12) 

We can now define a fictious rather then reconstructed velocity field v on the surface image in by taking the 
time derivative of each particle's positions given by x°. The vector field can also be recovered directly from qab — 
^\ A B) by the inverse of V^. We can then use v to evaluate the coarse-grained as tensors in E^ 

{Z{ab)) ^ [ V(^anb)da (13) 

J dDt 

(Vo denotes here the volume enclosed by dDt in E''). Because of ([8]) these integrals are independent of the vector field 
V chosen, but they do depend on the orientation of the embedded surface and under gauge transformations (jl2p they 
transform according to 

(Zab) R^a R'^^b {Zed)- 
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m 



FIG. 2: Canonical isometry Qp : TpT,t 



In order to get rid of that gauge dependence we need to map (Zab) back to the spacetime. Each isometric embedding 
ft defines canonical isometrics gp between the tangent spaces TpT,t to the constant time slices at points p £ dCt and 
E^. Any vector X eTpY.t can be decomposed into the tangent part to dCt and the projection to the outward pointing 
normal 

X' = X]. + Xnn\ 

We then take 

gp{XY = {ft.XTT +Xnm\ 

where m" is the outward pointing normal to the image of the surface at f{p) (see fig. [2). gp preserves the metric 
tensor and can be used to map geometric objects from TpSf to or back. In particular, one can pull {Zab) back to 
each point p in dCt by 5* 

{Z^,)p^g;({Z)) . 

It is straightforward to verify that the tensor field {Zij ) constructed on dCt this way is insensitive to transformations 
([T^ and thus completely independent of the choice of isometric embeddings /(. 

If the antisymmetric part of Z^^, i.e. the vorticity, does not vanish, it can be coarse-grained as well, although in 
a different way. Note that the 3+1 splitting cannot be orthogonal in that case and the four-velocity field u^^ has a 
non-vanishing projection to the constant time slice, denoted by w' . We define now the antisymmetric part of (Zab) 
as 

{Ziab]) = ^ I g{w)ianb]dcr. (14) 

^0 J dDt 

Note however that, in contrast to the rotating case, in a spacetime filled with rotational dust it is difficult to single 
out a preferred 3+1 splitting. Unfortunately, just like in the Buchert's formalism, the values of (Zij) depend on the 
splitting one chooses in general. Throughout the rest of the paper we will only deal with the non-rotating dust and 
in that case {Z^ij]) = 0. 

A. The coarse— graining procedure in the narrow tube limit and in FLRW metrics 

Outside the Newtonian setting we loose the equivalence of (Zij) to the volume average of the local velocity gradient. 
Nevertheless the newly defined quantities turn out to have the correct limiting behavior, namely if we shrink the coarse- 
graining domain to a single point, we recover the local Zij at that point. More precisely, consider a one-parameter 
family of tubes described by equation 




(15) 
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with a positive parameter A. As A ^> 0, the tube shrinks down to the single worldhne = 0. One can prove that 
for smaU A the values of (Zij), including the antisymmetric part if the vorticity doesn't vanish, can be expanded at 
every point p of the tube in powers of A as 

(Z,,)= Z,,|^.^0 + O(A). (16) 

Intuitively one can understand this fact in the following way: expressions (jl3p and ^ are perfectly equivalent if the 
background metric is flat and if the velocity gradient is constant over the coarse-graining domain. In general none of 
these conditions is satisfied on Sj, but the more the domain under consideration shrinks, the less it "feels" the back- 
ground curvature and the inhomogeneities in Zij. As it shrinks down to a point, both curvature and inhomogeneities 
drop out and (Zij) is equal to the volume average of the local Zij over a small domain around the point — 0. It 
is now straightforward to prove that such volume average approaches then the value of Zij at = 0. For a rigorous 
argument see Q. 

In order to test the coarse-graining procedure we may also apply it to known homogeneous cosmological solutions. 
For a dust-filled Friedman-Lemaitre-Robertson- Walker metric in the natural 3+1 splitting the it yields the expected 
answer for any admissible surface: 

{Zab) =H{t)Sab, (17) 

i.e. the coarse-grained velocity gradient consist only of the scalar part given by the Hubble parameter. This result 
is easy to see if we notice that the metric qab induced on the surface, like any other object in a FLRW solution 
comoving with the dust, undergoes a homogeneous rescaling in time described by — 211 (t) qab- The same metric 
variation is induced on the image of the embedded surface by the homothetic vector field w° = H{t) a;° in E^. It IS 
now straightforward to verify that the restriction of this vector field to dDt gives P7)) when plugged into ([T^ . This 
result holds for open, closed and fiat FLRW metrics alike. 



IV. TIME EVOLUTION EQUATION FOR THE COARSE-GRAINED QUANTITIES 

In order to derive the time evolution equations for the coarse-grained velocity gradient, one must first introduce a 
preferred time derivative of a tensor field in E^, which will serve as a generalization of the covariant derivative along 
the geodesic in (fTOl) . We define it via 



d_ 
di 



where 



Wab = - -TT / W[a " 6] dcT. 

^0 JdDt 

T> preserves the Euclidean metric in E^ and has the right narrow tube limit, i.e. it reduces to V„ along the limiting 
geodesic 0. 

We now fix a tube of particle worldlines Ct, given by functions £,'^{9^), whose sections dCt are admissible surfaces. 
We find the isometric embeddings and calculate {Zab)- This allows for decomposing the local Zij^ pushed forward to 
dDt C E'^ via ^p, and the velocity field into the the coarse-grained, large scale part and local inhomogeneities 

gp*{Z)ab = {Zab)+5Zab (18) 

= {Z\)x'' + 5v\ (19) 

both defined on dDt only. 

For the sake of brevity we introduce the following notation for surface integrals of the type of ^ and p^ : 

^[Xa]b^^ I XaUbAn. 

JdDt 

We will also denote by 7?. [ • ] the unique inverse of satisfying 
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for any tab, plus an irrelevant condition fixing the constant part for example fixing the surface's centroid at the 
origin. Now the definition of the coarse-grained Zij can be written down in a slightly more compact manner 



= N-TZ [qabU = A/"^ 2Z(,,) e^A ^\b 



(20) 



in which we have used the time derivative of ([TT|) to substitute '^Zuj-^ C a^'' b Qab- Note that the combination 



j\fn 

St. 



plays in our formalism the role of the coarse-graining operator for symmetric tensors of rank 2 on 
The evolution equation for (Zab) takes now the form of 

T>{Zab) = —{Zac) {Z^b) ^ {RaObo) + Bab + Bab, (21) 

where 

(RaObo) = J^Ti- '^RiOjO C A^^ B 

L ' ' J (ab) 

is obviously the coarse-grained contraction of the Riemann tensor with u'^, considered as a symmetric tensor on Et 
0. We can see that two new terms have appeared in comparison with (jlOp . The first one is the symmetrized version 
of the familiar Newtonian backreaction term ([5|) : 



Bab = M[Sv'^JV(a],^-J^ 



Sv' 



(22) 



The second one is entirely relativistic and has a more complicated structure 



B„ 



4(Z,,) (6Z''^^x\^A -^^\a) x\b) 



ab 



Nn [2 {5Z,, bZ\ x\a x\b - ^v\a Sv'^b Scd)]^, 



Just like in the Newtonian case, both backreaction terms are surface integrals divided by volume. In contrast to 
(|22p however, Bab involves linear terms in perturbations. Both Bab and Bab together measure the influence of 
inhomogeneities on the time evolution of the coarse-grained expansion and shear. 
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